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The growing recognition that entanglement is not exclusively a quantum property, and does not 
even originate with Schrodinger’s famous remark about it [Proc. Camb. Phil. Soc. 31 , 555 (1935)], 
prompts examination of its role in marking the quantum-classical boundary. We have done this by 
subjecting correlations of classical optical fields to new Bell-analysis experiments, and report here 
values of the Bell parameter greater than B = 2.54. This is many standard deviations outside the 
limit B = 2 established by the Clauser-Horne-Shimony-Holt (CHSH) Bell inequality [Phys. Rev. 
Lett. 23 , 880 (1969)], in agreement with our theoretical classical prediction, and not far from the 
Tsirelson limit B = 2.828.... These results cast a new light on the standard quantum-classical 
boundary description, and suggest a reinterpretation of it. 


Introduction: For many decades the term “entanglement” has been attached to the world of quantum mechanics 
[l]. However, it is true that non-quantum optical entanglement can exist (realized very early by Spreeuw Q) and 
its applications have concrete consequences. These are based on entanglements between two, or more than two, 
degrees of freedom, which are easily avalable classically Multi-entanglements of the same kind are also being 

explored quantum mechanically Q. Applications in the classical domain have included, for example, resolution of 
a long-standing issue concerning Mueller matrices Q, an alternative interpretation of the degree of polarization Q, 
introduction of the Bell measure as a new index of coherence in optics 0 , and innovations in polarization metrology 
11]. Here we present theoretical and experimental results extending these results by showing that probabilistic classical 
optical fields can exhibit violations of the Clauser-Horne-Shimony-Holt (CHSH) Bell inequality 12] of quantum 
strength. This is evidence of a new kind that asks for reconsideration of the common understanding that Bell 
violation signals quantum physics. We emphasize that our discussion focuses on non-quantum entanglement of non- 
deterministic classical optical fields, and does not engage issues such as non-locality that are important for some 
applications in quantum information. 

The observations and applications of non-quantum wave entanglement noted above @ 11§0 exploited non- 
separable correlations among two or more modes or degrees of freedom (DOF) of optical wave fields. Nonseparable 
correlations among modes are an example of entanglement fl3j . but are not enough for our present purpose. In 
addition, we want to conform to three criteria that Shimony has identified for Bell tests [I], facts of quantum Nature 
that must be satisfied when examining possible tests of the quantum-classical border. Fortuitously, the ergodic 
stochastic optical fields of the classical theory of partial coherence and partial polarization (see Wolf [l5() satisfy these 
criteria fully (see Suppl. Materials 16]), and we have used such fields as our test bed. 

Background Theory: We will deal here only with the simplest suitable example, the theory of completely unpolar¬ 
ized classical light, and have explained elsewhere (see 0. a ) the generalizations needed to treat partially polarized 
fields, which lead to the same conclusions. In all cases there are only two degrees of freedom (DOF) to deal with, 
namely the direction of polarization and the temporal amplitude of the optical field. In both classical and quantum 
theories these are fundamentally independent attributes. An electric field, for a beam travelling in the z direction, is 
written 


E(t) = xE x (t) + yE v (t). 


( 1 ) 


In the classical theory of unpolarized light [18] the optical field’s two amplitudes E x and E y are statistically completely 
uncorrelated, and are treated as vectors in a stochastic function space. A scalar product of the vectors in this 
space corresponds physically to observable correlation functions such as (E x E y ). For unpolarized light we have 
( E X E X ) = ( E y E y ), and ( E x E y ) = 0. 

Now it is possible to talk of entanglement of the classical field. This is because entangled states are superpositions 
of products of vectors from different vector spaces, whenever the superpositions can’t be rearranged into a single 
product that separates the two spaces [I] . Looking again at ©, we see that this is the case because we’ve taken E to 
be unpolarized. That is, by the definition of unpolarized light, there is no direction u of polarization that captures 
the total intensity, so E{t) can’t, for any direction u, be written in the form E(t) = uF(t ), which would factorably 
separate the polarization and amplitude DOF 0 - 

Beyond its probabilistic indeterminacy, the E in © has other quantum-like attributes - it has the same form as a 
quantum state superposition and can be called a pure state in the same sense, more precisely a two-party state living 
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in two vector spaces at once, polarization space for x and y , and infinitely continuous stochastic function space for 
E x and E y . 

The Bell inequality most commonly used for correlation tests is due to Clauser, Horne, Shimony and Holt (CHSH) 
fl3 |. It deals with correlations between two different DOF when each is two-dimensional. The Schmidt Theorem of 
analytic function theory [20| ensures two-dimensionality, by guaranteeing that among the infinitely many dimensions 
available to the amplitudes in ©, only two dimensions are active. This is a consequence arising just from the fact 
that the partner polarization vectors x and y live in a two-dimensional space. 

For convenience, we introduce e, the field normalized to the intensity I = (E X E X + E y E y ): 

e(t) = E{t)/Vl = {xe x {t) + ye y (t)}, ( 2 ) 


where now (e • e) = (e x e x + e y e y ) = 1. 

For some simplification in writing, we will use Dirac notation for the vectors without, of course, imparting any 
quantum character to the fields. The unit polarization vectors x and y will be renamed as x —> |zti) and y — > |zi 2 ) 
and the unit amplitudes will be rewritten e x —> |/i) and e y —> I/ 2 }. If desired, the Dirac notation can be discarded 
at any point and the vector signs and hats re-installed. For the case of unpolarized light we have (wi|rt 2 ) = 0 and 
(/ 1 I/ 2 ) = 0. Unit projectors in the two spaces take the form 1 = |zzi)(zzi| + \u 2 ){u 2 \ and 1 = |/i)(/i| + |/ 2 )(/ 2 |- In 
this notation, and in the original notation for comparison, the field takes the form 

E/Vi = x e x + y e y 

= l e ) = (l'»i)l/i) + M|/2 ))/v / 2. (3) 

In this notation the field actually looks like what it is, a two-party superposition of products in independent vector 
spaces, i.e., an entangled two-party state (actually a Bell state). Here the two parties are the independent polarization 
and amplitude DOF. 

The notation for a CHSH correlation coefficient C(a,b) implies that arbitrary rotations of the unit vectors \uj) 
and | fk){j, k = 1, 2) through angles a and b can be managed independently in the two spaces. An arbitrary rotation 
through angle a of the polarization vectors |zti) and \u 2 ) takes the form 


| zz“) = cosa|wi) — sina|zi 2 ) and 
IZZ 2 } = sinajui) + cosa|zt 2 ). 

For function space rotations we have | f\) and |/|) defined similarly: 

\fi) = cos b\fi) - sin b\f 2 ) and 
I/I) = sin fo| / 1 ) +COS&I/ 2 }, 


(4) 

(5) 


where the rotation angles a and b are unrelated. 

Next, the correlation between the polarization (zz) and function (/) degrees of freedom is given by the standard 
average 


C{a,b) = {e\Z u (a)®Z f (b)\e), (6) 

where Z is shorthand for the difference projection: Z u (a) = |zz“)(zz“| — Izt^Xzz^l, analogous to a a z spin operation. 
C(a, b) is thus a combination of four joint projections such as: 

Pi rM) = ( e l (l u i)l/i)(/i l( u il) l e ) (7) 

= (/il( u il e ) 2 ■ 


This is all classical and all of the correlation projections Pjk{a, b) with j,k = 1,2 have familiar roles in classical optical 
polarization theory [l8( . 

Gisin Q observed that any quantum state entangled in the same way as the classical pure state © will lead to 
violation of the CHSH inequality, which takes the form B < 2, where 


B = | C(a, b ) - C(a', b) + C(a, b') + C{a ', b')\. 


(8) 
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The same result will be found here, as one uses only DOF independence and properties of positive functions and 
normed vectors to arrive at it (see details in Suppl. Matls. [lgj]). We note again that the issue of entanglement itself 
is pertinent to the discussion, but the usefulness of entanglement as a resource for particular applications is not. Thus 
we have reached the main goal of our theoretical background sketch. This was to demonstrate the existence of a 
purely classical field theory that can exhibit a violation of the CHSH Bell inequality. 

Experimental Testing: The remaining task is to show that experimental observation confirms this theoretical 
prediction, in effect shifting one’s interpretation of tests of the quantum-classical border by showing that, along with 
quantum fields, classical fields conforming to the Shimony Bell-test criteria are capable of Bell violation. In order to 
make such a demonstration, a classical field source must be used. This means one producing a field that is quantum 
mechanical (since we believe all light fileds are intrinsically quantum), but a field whose quantum statistics are not 
distinguishable from classical statistics. This is only necessary up to second order in the field because the CHSH 
procedure engages no higher order statistics. Such sources are easily available. Since the earliest testing of laser light 
it has been known that a laser operated below threshold has statistical character not distinguishable from classical 
thermal statistics. So in our experiments we have used a broadband laser diode operated below threshold. 

Our experiment repeatedly records the correlation function C(a,b) defined in © for four different angles in order 
to construct the value of the Bell parameter B. This is done through measurements of the joint projections Pjk(a, b). 
We will describe explicitly only the recording of Pu(a, b), identified in (|7|). but the others are done similarly in an 
obvious way. In the classical context that we are examining, the optical field is macroscopic and correlation detection 
is essentially calorimetric (i.e., not requiring or employing individual photon recognition). 

Polarization Tomography: The first step is to determine tomographically the polarization state of the test field. 
A polarization tomography setup is shown in Fig. 1. Using a polarizing beam splitter and half and quarter wave 
plates to project onto circular and diagonal bases, the Stokes parameters (Si, S 2 , S 3 ), relative to So = 1, are found 
to be (—0.0827, —0.0920, —0.0158), providing a small non-zero degree of polarization equal to 0.125. This departure 
from zero requires a slight modification of the theory presented above (see the Supplemental Material) and reduces 
the maximum possible value of B able to be achieved for our specific experimental field to B = 2.817, below but close 
to B = 2\[2 = 2.828..., the theoretical maximum for completely unpolarized light. 

Experimental Bell Test: The experimental test has two major components, as shown in Fig. [TJ a source of the 
light to be measured, and a Mach-Zehnder (MZ) interferometer. The source utilizes a 780 nrn laser diode, operated in 
the multi-mode region below threshold, giving it a short coherence length on the order of 1 mm. The beam is assumed 
statistically ergodic, stable and stationary, as commonly delivered from such a multimode below-threshold diode. It is 
incident on a 50:50 beam splitter and recombined on a polarizing beam splitter after adequate delay so that the light 
to be studied is an incoherent mix of horizontal and vertical polarizations before being sent to the measurement area 
via a single mode fiber. A half wave plate in one arm controls the relative power, and thus the degree of polarization 
(DOP). Quarter and half wave plates help correct for polarization changes introduced by the fiber. 

In Fig. Q] the partially polarized beam entering the MZ is separated by a 50:50 beam splitter into primary test 
beam |E) and auxiliary beam |E). The two beams inherit the same statistical properties from their mother beam 
and thus both can be expressed as in Eq. ©, with intensities I and I. The phase of the auxiliary beam |E) is shifted 
by an unimportant factor i at the beam splitter. 

To determine the joint projection Pn(a, b) of the test beam |E), the first step is to project the field to obtain 
|E“ ) = |-u“)(u“|E). This can be realized by the polarizer labelled a on the bottom arm of the MZ. The transmitted 
beam retains both |/) components in function space: 

\K) = V^\<)(cii\f b 1 )+c 12 \f b 2 )), (9) 

where /“ is the intensity, and cn and C 12 are normalized amplitude coefficients with |cn| 2 + |ci 2 1 2 = 1. Here cn 
relates to Pu in an obvious way: Pn(a,b) = /“|cn| 2 //. One sees that the intensities I and /“ can be measured 
directly, but not the coefficient C\\. 

For Pn(a,b) our aim is to produce a field that combines a projection onto | f\) in function-space with the |u“) 
projection in polarization space. The challenge of overcoming the lack of polarizers for projection of a non-deterministic 
field onto an arbitrary direction in its independent infinite-dimensional function space is managed by a “stripping” 
technique p3] applied to the auxiliary E field in the left arm. We pass E through a polarizer rotated from the initial 
|iti) — | U 2 ) basis by a specially chosen angle s, so that the statistical component I/I) is stripped off. The transmitted 
beam |Ef) then has only the |/|) component, as desired: |Ef) = i\/T? |uf)| f, ). Here 7f is the corresponding intensity 
and the special stripping angle s is given by tans = («i/K 2 )tan& (see [3L ll6l] 1. 

The function-space-oriented beam |Ef) is then sent through another polarizer a to become |E“) = |w,J)( u i |E“) = 
*y / /j*j'u“)l/i); where /“ is the corresponding intensity. Finally, the beams |E“) and |E“) are combined by a 50:50 
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FIG. 1: The experimental setup consists of a source of unpolarized light and a measurement using a modified Mach-Zehnder interferometer. 
Half and quarter wave plates (HWP, QWP) control the polarization of the source. All beam splitters are 50:50 unless marked as a polarizing 
beam splitter (PBS). Intensities needed for obtaining the required joint projections are measured at detector Dl. Shutters S independently 
block arms of the interferometer in order to measure light through the arms separately. A removable mirror (RM) directs the light to 
a polarization tomography setup, where the orthogonal components of the polarization in the basis determined by the wave plate are 
measured at detectors D2 and D3. 


beam splitter which yields the outcome beam |Ef) = (|E“) +i|E“))/\/2- The total intensity if of this outcome beam 
can be easily expressed in terms of the needed coefficient cn. 

Some simple arithmetic will immediately provide the joint projection Pn(a, 6) in terms of various measurable 
intensities: 

Pn(a, b) = (2 If -If- If) 2 /4llf. (10) 

Other Pjk{a,b) values can be obtained similarly by rotations of polarizers a and s. To make our measurements, 
polarizers a were simultaneously rotated using motorized mounts, while the third polarizer s was fixed at different 
values in a sequence of runs. 

Results: For each angle, measurements were made at detector Dl for the total intensity I T , and the separate 
intensities from each arm I a and I a by closing the shutters S alternately. In this way, the measurements of the 
polarization space and statistical amplitude space are carried out separately. From these measurements the needed 
correlations C(a, b) were determined and Eq. © used to evaluate the CHSH parameter B. 

Fig. [5]shows C(a, b) obtained by measuring the joint projections Pjk{a, b) for a complete rotation of polarizer a, with 
different curves corresponding to b (and thus s ) fixed at different values. It is apparent from the near-identity of the 
curves that, to good approximation, the correlations are a function of the difference in angles, i.e. C(a, b ) = C(a — b). 
The maximum value for B can then be found straightforwardly from any one of the curves in Fig. [2j Among them 
the smallest and largest values of B (obtained for curves 1 and 4), are 2.548 ± 0.004 and 2.679 ± 0.007. 

To be careful, we note that in our experiments the field was almost but not quite completely unpolarized, thus not 
quite the same field sketched in the Background Theory paragraphs. Thus we couldn’t expect to get the maximum 
quantum result B = 2yp2 = 2.828... for the Bell parameter, but the values achieved also present a strong violation. 
The background theory is mildly more complicated for partially polarized rather than unpolarized light, but when 
worked out for the degree of polarization of our light beams (see [3[ and the Suppl. Matl. |16j) it supports the values 
we observed. 
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polarizer angle a 


FIG. 2: Plots of the correlation functions C(a,b) obtained by rotating polarizer a in the polarization space and holding angle b in the 
function space constant. Curves 1 to 4 correspond to different fixed values of b separated by 7t/4. The invariant cosine function required 
to violate the Bell inequality is clearly present. Error bars are included but scarcely visible. 


Summary: In summary, we first sketched the purely classical theory of optical beam fields (JTJ) that satisfy the 
Bell-test criteria of Shimony [lj, [l6[ . Their bipartite pure state form shows the entanglement of their two independent 
degrees of freedom [22j. The classical theory defines them as dynamically probabilistic fields, meaning that individual 
field measurements yield values that cannot be predicted except in an average sense, which is another feature shared 
with quantum systems but also associated for more than 50 years with the well-understood and well-tested optical 
theory of partial coherence 18]. Our theoretical sketch for the simplest case, unpolarized light, indicated that such 


fields or states are predicted to possess a range of correlation strengths equal to that of two-party quantum systems, 
that is, outside the bound B < 2 of the CHSH Bell inequality and potentially as great as B = 2\j2. In our experimental 
test we used light whose statistical behavior (field second-order statistics) is indistinguishable from classical, viz., the 
light from a broadband laser diode operating below threshold. Our detections of whole-beam intensity are free of the 
heralding requirements familiar in paired-photon CHSH experiments. Repeated tests confirmed that such a field can 
strongly violate the CHSH Bell inequality and can attain Bell-violating levels of correlation similar to those found in 
tests of maximally entangled quantum systems. 

One naturally asks, how are these results possible? We know that a field with classically random statistics is a 
local real field, and we also know that Bell inequalities prevent local physics from containing correlations as strong as 
what quantum states provide. But the experimental results directly contradict this. The resolution of the apparent 
contradiction is not complicated but does mandate a shift in the conventional understanding of the role of Bell 
inequalities, particularly as markers of a classical-quantum border. Bell himself came close to addressing this point. 
He pointed out Q that even adding classical indeterminism still wouldn’t be enough for any type of hidden variable 
system to overcome the restriction imposed by his inequalities. This is correct as far as it goes but fails to engage the 
point that local fields can be statistically classical and exhibit entanglement at the same time. For the fields under 
study, the entanglement is a strong correlation that is intrinsically present between the amplitude and polarization 
degrees of freedom, and it is embedded in the field from the start (as it also is embedded ab initio in any quantum 
states that violate a Bell inequality). The possibility of such pre-existing structural correlation is bypassed in a 
CHSH derivation. Thus one sees that Bell violation has less to do with quantum theory than previously thought, but 
everything to do with entanglement. 
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SUPPLEMENTARY MATERIAL 

Shimony’s Bell-Test Conditions: In his extended analyses jl] of Bell inequalities and their testing, Shimony 
recognized that in order to deserve serious attention, an alternative non-quantum theory entering what is considered 
a quantum domain (in our case, the domain of Bell test violation) needs to embrace in some way aspects of Nature 
that appear completely random, i.e., purely probabilistic. These aspects are dealt with by quantum theory in well 
known ways, and this is the reason Bell once raised the issue of classical indeterminism Q, but without going as 
far as classical entanglement. Shimony summarized these considerations by naming three key features, all of which 
should be considered, in his words, “... as established parts of physical theory: (I) In any state of a physical system 
S there are some eventualities which have indefinite truth values. (II) If an operation is performed which forces an 
eventuality with indefinite truth value to achieve definiteness ... the outcome is a matter of chance. (Ill) There are 
‘entangled systems’ (in Schrodinger’s phrase) which have the property that they constitute a composite system in a 
pure state, while neither of them separately is in a pure state.” Here by eventualities Shimony means measurement 
outcomes. We have identified the electric field as it is dealt with in the standard classical theory of partial coherence 
and partial polarization as a physical system satisfying all three conditions. 

Partially Polarized Fields: Classical statistical light of any degree of polarization can be treated exactly as the 
unpolarized light in the text, with a small change. One needs to insert the parameters that allow the two orthogonal 
components provided by the Schmidt decomposition to have different intensities. These parameters, rei and « 2 , are 
the (real) Schmidt eigenvalues. They satisfy k\ + re| = 1, and both equal l/\/2 in the completely unpolarized case 
treated in the text. Any intensity-normalized partially polarized field can then be written as in text Eqn. (3), but 
with m and «2 attached Q: 


|e) = Kl|ztl)|/l) +K 2 \u 2 )\f 2 )- (11) 

It is clear that the field is entangled between polarization and amplitude unless one of the re’s is zero, in which case |e) 
is plainly separable. All of the formulas for unpolarized light will change, but only to the extent that the presence of 
the re’s requires. The conventional degree of polarization P is given in terms of the re’s by P = |ref — re|| |^|, and this 
can be used to find the Schmidt coefficients. For the experimental classical statistical optical light field we obtained 
Ki, K 2 = 0.750, 0.661. Another interesting formula is the full result for the partially polarized Bell parameter: 

B = cos2a(cos26 — cos 2 &') 

+ cos 2 a' (cos 2 b + cos 2 b') 

+ 2 reire 2 {sin 2 a(sin 26 — sin 21 /) 

+ sin 2 a' (sin 26 + sin 2 b')}. ( 12 ) 

One sees two uncomplicated limits: when either re is zero (no entanglement) no result higher than 3 = 2 can be 
achieved, and when rei = re -2 = 1/V2 (maximal entanglement) the Tsirelson bound can be reached: B = 2\/2. Also, 
if one follows Gisin’s approach :5| by choosing the rotation angles as a = 0, a' = 7 r/ 4 , and cos2b = —cos2b' = 
l/i/l + 4 |reiK 2 | 2 , the Bell parameter becomes 


B = 2/ y/l + 4|kik 2 | 2 . (13) 

Apparently, there will be a Bell violation {Bl2) as long as reiK 2 ^ 0, i.e., when there is non-zero entanglement. 

Basis stripping, rotation, and projection in amplitude function space: To observe correlation, it is essential 
to be able to access and measure both polarization and amplitude function degrees of freedom. Unfortunately, unlike 
the polarization degree of freedom, there is no systematic technology working directly in the infinite-dimensional 
amplitude function space that can project a non-deterministic field onto an arbitrary basis |/f) in that space. This 
requires an innovation using an indirect measurement such as incorporated in the experimental setup sketched in Fig. 
1 of the text. In our setup we employ an auxiliary beam that contains only the |/f) basis to interfere with the primary 
test beam (which in general contains both |/*) and |/|)). By interference measurements one is able to obtain the |/*) 
information from the test beam, and finally determine the needed Cn value given in equation (9) of the main text. 

This section shows specifically how to “strip” off and rotate a basis in the statistical amplitude function space of 
the auxiliary field |E), which shares the same statistical properties as the primary test beam. For generality, we take 
the auxiliary beam as initially in the form of Eq. m with arbitrary rei, K2. Such a beam can always be rewritten in 




the rotated amplitude function space basis |/f), |/I), i.e, 


|E) 



(ki cos b\u\) 


k 2 sin b\u 2 )) |/i) 


+ (ki sin b\u\) + n 2 cos6|u 2 }) \f 2 ) 


(14) 


One notes from the second term of the equation that a properly chosen polarizer that blocks completely the 
polarization-space component Kisin&|rti) + n 2 cos b\u 2 ) will effectively strip off the amplitude function space basis 
vector |/|). 

Such a stripping polarizer |uf) is oriented with respect to the polarization space basis |ui), \u 2 ) by a specific angle 
we have called s, i.e., 


\u{) = coss|ui) — sins|u 2 ). (15) 

Then the stripping condition is directly given as 

(u^l («i sin5|rti) + k 2 cos6|m 2 }) = 0, (16) 

which specifies the rotation angle s according to 

tan s = (ki/k 2 ) tan b, (17) 

so s is determined by the values of n\ and n 2 for any choice of rotation angle b in the amplitude function space. 

As a result of this stripping polarizer, the beam (TT4l) becomes 

|Ef) = K)(«f|E) 

= \//(ki cos&coss + k 2 sin5sins)|it*}|/ i), (18) 

where the amplitude function space component |/|) of the transmitted beam is completely stripped off. 

One notes immediately that arbitrary rotations of the function space basis |/f) (i.e., arbitrary change of angle b) for 
the auxiliary beam can be effectively realized by the change of angle s through Eq. m Consequently, as described 
in the main text, by further interference of this auxiliary beam with the primary test beam, one is able to project the 
test beam onto any function space basis |/j ) ) with an arbitrary angle b. This is exactly how we make measurements 
in the amplitude function space of the statistical light beam. 
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